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Abstract 

For a machine to be useful in practice, it preferably has to meet two requirements: namely, 
(i) to be able to perform work under a load and (ii) its operational regime should ideally not 
depend on the time at which the machine is switched-on. We devise a minimal setup, consisting 
of two atoms only, for an ac-driven quantum motor which fulfills both these conditions. Explicitly, 
the motor consists of two different interacting atoms placed into a ring-shaped periodic optical 
potential - an optical "bracelet" -, resulting from the interference of two counter-propagating 
Laguerre-Gauss laser beams. This bracelet is additionally threaded by a pulsating magnetic flux. 
While the first atom plays a role of a quantum "carrier" , the second serves as a quantum "starter" , 
which sets off the "carrier" into a steady rotational motion. For fixed zero-momentum initial 
conditions the asymptotic carrier velocity saturates to a unique, nonzero value which becomes 
increasingly independent on the starting time with increasing "bracelet" -size. We identify the 
quantum mechanisms of rectification and demonstrate that our quantum motor is able to perform 
useful work. 
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I. INTRODUCTION 



The development of laser cooling techniques and the creation of the Bose-Einstein con- 
densates (BEC) present landmark successes of modern experimental physics^. Since 1995-, 
ultracold atoms rapidly become a popular toolbox to explore the quantum world. Subse- 
quent experimental studies can be (conditionally) divided into three stages. The first one 
was aimed to model the electronic behavior in solids with ultracold atoms trapped in peri- 
odic optical potentials. These potentials, so called optical lattices thus implement the idea of 
"quantum simulators", proposed by Feynman in 1982- . Incarnation of paradigmatic quan- 
tum models (Bose- and Fermi-Hubbard systems, Tonk-Girardeau gas and etc)^ and their 
intensive experimental studies mark the second stage. The recent advent of a third wave, 
namely "to make use of it", triggered a search for applications which led to single atom 
qubits^, atom optical clocks-, cold atom interferometry-l, as well as cold atom gyroscopes^. 

An electric motor, i.e. a device that converts electrical energy into mechanical work by 
setting a working body into a linear or rotational motion, is a archetype example of a useful 
physical system. For nearly two centuries, since the invention of the first electrical motor-, 
a never-ending continuous miniaturization has already passed the microscale leveU^ and 
entered the nano-scale worlci^. Yet, this process alone does not parallel the transition from 
the classical to the quantum world: the operational descriptions of all existing electrical, 
nano-sized motors rest on classical concepts^^. 

With this present work, we provide a detailed analysis for an electric quantum motor 
proposed recently in Ref^^. This motor is made of two ultracold atoms only that are trapped 
in a deep ring-shaped one-dimensional optical lattice, - an optical "bracelet". The blueprint 
for such an underlying trapping potential was proposed recently^^ and a first experimental 
realization has been reported in Ref^^. We employ this setup to devise an engine which 
works as a genuine ac-driven quantum motor. We identify the quantum mechanisms which 
yield the modus operandi for a motor device. Moreover, we discuss parameter values suitable 
for the realization of the quantum engine with present-day experimental setups. 
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FIG. 1: Sketch of the quantum machine: Two ultracold neutral atoms are loaded into an optical 
"bracelet" - a ring-shaped optical lattice. Both atoms interact locally with each other, while only 
one carrier (the one with an arrow) is magnetically driven. 

II. THE MOTOR SETUP 

Figure 1 is a sketch of our quantum machine device. The optical potential, which results 
either from the interference of a Laguerre- Gauss (LG) laser beam with a plane wave^- or, 
alternatively, of two collinear LG beams with different frequencies^' is capable of trapping 
two interacting atoms. The first atom, termed "carrier", c, is assumed to be charged and 
is driven by the time-dependent magnetic flux piercing the bracelet^. The second atom, 
termed "starter", s, is neutral thus remains unaffected by the magnetic flux, but interacts 
locally, e.g. by means of atom-s-wave scattering, with the carrier when both atoms share 
the same site of the optical lattice. 

We next assume that both atoms are loaded into the lowest energy band of a deep, ring- 
shaped optical potential with L lattice sites with the lattice constant d. The time-dependent 
homogeneous vector potential A{t) does not induce any appreciable transitions between the 
ground band and the far separated excited band(s). 

The total Hamiltonian of the system reads 



tot — 



int ) 



(1) 



where the time-dependent Hamiltonian for the carrier is given by 
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and the starter Hamiltonian Hg reads 

Hs = -^ ^5^|/. + l)(/.|+H.c.j (3) 

Here, Jc and Jg are the corresponding hopping strengths to neigbouring sites, which depend 
on the atom mass, M, and the potential depth, Vq. Here, for both atoms, we have assumed 
the hmit of a deep periodic potential, also refereed to as the tight-binding model. This limit 
is well justified for a potential amplitudes of Vq ^ 5Eq, where the higher-order tunneling 
amplitude is less than 10% of the leading first-order (i.e., the next neighbour) tunneling 
amplitude^. Eq = h'^TT^/2Md^ is the "recoil" energy. The salient carrier-starter on-site 
interaction is 

L 

Hint = W J2 ^lc,ls\Qm^\ls){ls\, (4) 

where W denotes the interaction strength. Throughout the remaining, we use periodic 
boundary conditions; i.e., |L -|- 1) = The dimension of the complete Hilbert space is 
M = L?' . The total system is described by wave function \%l){t)). The scale of the motor 
current will be measured in units of the maximal group velocity i^o = Jcd/h. 

To conclude with the setup, we specify that the driving is switched on at the time t^, so 
that the vector potential has the form 

A{t-to)=x{t-to)A{t). (5) 

where xi^ — ^o) is the step function, and A{t) is defined on the entire time axis, t e 
(-00, +oo). 



III. DC-QUANTUM CURRENT 

The mean carrier current is given as the speed of the motor by means of the velocity 
operator: Vc(t;to) = i/h[Htot(t),Xc]- With = d'^ihlh) {h\ ® Is, one finds Vcit'^io) = 
—i{vo/2) (^Ylfc=i e*"^^*'*"'' l^c + 1) {lc\ — H.c.) ® !«. In the quasimomentum representation with 
= X]n=i 6xp(min)|r;,), its quantum expectation Vc{t]to) = {ip{t)\vc{t]to)\ijj{t)) reads 

L 

Vc{t; to) =Vo'^ p^^ (t; to) sin (^ki + A{t; to)) , (6) 
1=1 
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wherein ki = 2tiI /L is the single particle quasimomentum and where we indicated explicitly 
the parametric dependence on the start time to. Further, Pf^it^to) = \{4'{'t)\'^h^ '^iJl'^j 
where \ki^,k,iJ = \k,i^) ® is the quasimomentum distribution for the carrier. The 

asymptotic steady state regime of the motor can be characterized by the dc-component of 
the averaged velocity; i.e., 

1 /•* 

Vc{to) := hmt^oo- — — / Vcis;to)ds. (7) 
t - to Jto 

Without interaction between the particles, i.e., if W = 0, for an initially localized carrier 
with zero velocity not even a transient directed current does emerge^. This fact holds for 
any shape of the potential A{t). This situation mimics the one taking place in a single-phase, 
classical ac-motor: a periodically pulsating magnetic field would fail to put a rotor from rest 
into rotation, unless one applies an initial push via a starter mechanismi^. Similarly, a single 
quantum particle when initially localized on a single potential minimum of a periodically 
modulated ring-shaped potential doesn't acquire any momentum under an unbiased periodic 
driving force-^, while a constant bias, Asit) = ujst, induces Bloch oscillations only^. In our 
setup, it is the ingredient of the interaction with the second particle that takes over the role 
of a quantum starter. 

Notably, even for nonvanishing interaction, i.e. \W\ > 0, it is still not obvious how to set 
motor into rotation. Yet again, the seemingly "evident" solution - to apply a constant bias 
to the carrier - cannot solve the task. Due to the finite coupling with the starter system 
in Eq. ([3]), the corresponding vector potential, = UBt, may induce irregular Bloch 

oscillations with a resulting zero drift velocity only^^. In distinct contrast, we use here an 
unbiased time-dependent vector potential possessing a zero dc-component, A{T)dr = 
and being periodic in time, A{t + T) = A{t). 

Like for nanomechanical devices^^, the symmetry principles are of salient importance 
also for quantum engines: for a zero-momentum initial condition, an ac-force with time- 
reversal symmetry would launch the system - with equal probabilities - into a clockwise 
(rightward motion) or a counterclockwise rotation (leftward motion)^. Taking into account 
the quantum nature of the engine, the rotor will just spread symmetrically in both direction. 
Thus, the modus operandi requires a symmetry-breaking driving field, realized here with the 
harmonic mixing signal: 

A{t) = Aisin{ujt) + A2sin{2ujt + Q). (8) 
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where denotes the symmetry-breaking phase shift. The input ([8]) knowingly may induce 
a non-vanishing nonhnear response, the so-called ratchet effect^^^^^i^^^. 

IV. FLOQUET STATES AS TRANSPORT STATES 

The dynamics at times t > to of the time-periodic Hamiltonian ([1]) can be analyzed 
by using the Floquet formalism--. The solution of the eigenproblem: U{t,to)\(j)n{t', k)) = 
exp (-fe„(t - to)) \<Pn{t; k)), with the propagator U{t,to) = Texp (^-^ i^tot(r)j dr (T 
denotes the time ordering), provides the set of time-periodic Floquet states, with T = 2tt/u 
being the driving period, |0„(t + T; k)) = |0n(t; k)). Here, k = 27il/L with I = 1, . . . , L is 
the total quasimomentum of the Floquet state. Due to the discrete translation invariance of 
the system, the total quasimomentum is conserved during the time evolution, thus serving 
quantum number. 

In the absence of the driving, A(t) = 0, the motor setup ([UE]) possesses the continuous 
translational symmetry in time. In this case, the expansion coefficients of the initial wave- 
function ipito) in the system eigenbasis knowingly do not depend on time. On the contrary, 
eigenstates of a periodically driven system - the Floquet states - evolve in time, being locked 
by the external ac-field. Thus, the expansion of an initial wave-function over the Floquet 
eigenbasis depends on the start time to (E]), which determines the phase of the driving ac- 
field^^, i.e., |V'(^o)) = Zl^i Cn(to)|0„(to; /c)), with c„(to) = (0„(to; fc)|^/'(to)). Substitution of 
the above decomposition into ([7]) yields the result 

'^cito) = V't^n |c„(to)|^ ,Vn = — / Vn{t)dt . (9) 

n=l ^ -^to 

Here, Vn denotes the velocity expectation value of the n-th Floquet state ([6]). Because the 
Floquet states are periodic in time functions, the velocities U„ do not depend on to, and the 
dependence of the generated dc-current on the to solely stems from the coefficients c„(to). 
Since the system evolution is fully quantum coherent; i.e. there is no memory erasing induced 
by an environment, - the asymptotic current maintains the memory of the initial condition 
as encoded in the coefficients c„(to)^. 
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V. INPUT/OUTPUT CHARACTERISTICS 



The question is now, how can we control the motor? To answer this question, we used 
the symmetry analysis^S which allows us to predict an appearance of a certain dc current. 
Combining time-reversal operation and the complex conjugation applied to ([T]) with A{t) 
in the form of, one can prove the (anti-) symmetric dependence of U„ on for the Floquet 
states with A; = 0: (i) Un(— ©) = —Vn{Q), and (ii) U„(7r — 0) = U„(6). The first relation 
implies tJ„(0) = — U„(0) = 0. Then, for 6 = 0, the second relation gives tJn(7r) = 0. Thus 
the Floquet states with k = possess zero mean velocities at G = 0, vr. Furthermore, using 
a similar reasoning, one finds that the set of Floquet states with nonzero k can be ordered 
by the parity relation, which links eigenstates with opposite quasimomenta, 0„(t; —k; 9) = 
4>m(T — t] k; —9), yielding !;„(— 9) = —Vm{Q)- This implies that for a symmetric (in k) 
initial state and 9 = 0, vr, the contributions to the dc-current of Floquet states with opposite 
quasimomenta eliminate each other. The same holds true for a monochromatic driving ([8]), 
with A2 = 0^°. Shifting 9 away from 0, ±77 causes the decisive symmetry breaking and leads 
to the de-symmetrization of the Floquet states with k = and consequently will violate the 
parity between states with opposite signs of k. 

The emergence of a non- vanishing dc current at 9 7^ 0, vr induced by the coupling to the 
starter can be illustrated by the desymmetrization of velocities, pi^{t]tQ)sin + y4(t; tg) j 
(as a function of quasimomentum k), whose sum according to ([6]) yields the mean velocity 
of the individual Floquet states. For those with = 0, the desymmetrization happens more 
drastically since they do not produce any dc current at the symmetry point 9 = 0, ivr. 
On Fig. [21 we depict the instantaneous velocities, Pnit] to)sin (^k + A(t; to)j, setting to = 0, 
for a Floquet state with = 0, for 9 = 0, see Fig. [2]^a), and 9 = 7r/2, see Fig. [2](b), 
together with the resulting mean velocities Vn{t) Fig. [2]^c). The contributions from different 
quasimomenta for 9 = eliminate each other and give a periodically varying, but cycle- 
averaged zero current, see in Fig. EJj^c), solid line, while the desymmetrization for 9 7^ (b) 
results in non- vanishing dc current, cf. Fig. EJ^c), dashed line. 

The motor speed depends on the initial conditions, which define the contributions of 
different Floquet states to the carrier velocity (Q. We restrict our analysis to the initial 
state Itpito)) = L'^^'^\lc) ® h = 1, ■■■,L, in the form of the localized carrier (at Ic) 

and the uniformly "smeared", delocalized starter. Both particles assume zero velocities at 
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FIG. 2: (color online) The colormaps render velocities, Pnit', to) sin + A(t)j , whose sum accord- 
ing to ([6]) gives the mean velocity of an individual Floquet state with the total quasimomenta k = 
for (a) ^ = and (b) 6 = tt/2, as a function of quasimomentum k. (c) Velocity of the Floquet 
states, Vn{t), for ^ = (solid line) and 6 = ir/l (dashed line). Note that while v{t) for = 
oscillates and is non-zero except few points, its average over period, Vn, is strictly zero due to the 
symmetry Vn{t) = —Vn{T — t). The other parameters are = 0.1 x J^., = 0.5, A2 = 0.25, 
W = 0.2Jc, J^ = J^ = J, L = 16, and to = 0. 

t = to. The asymptotic velocity may exhibit a strong dependence on tcr^. We first discuss 
the results obtained after averaging over to, thus assigning a unique motor velocity value. 



for fixed system parameters. 

Figure [3] depicts the dependence of the average motor velocity on G. The results ob- 
tained by direct time propagation of the initial state and averaged over to (dashed line) are 
superimposed by those calculated via the Floquet formalism (solid line). Both curves 
show the expected symmetry properties Vc{Q) = —Vc{Q + n) = —Vc{—Q). The agreement 
between the two curves is satisfactory, although not perfect: This is so because the sharp 
peaks on the asymptotic motor velocity ([9]) and is due to the finite evolution time used in 
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FIG. 3: (color online) Averaged motor velocity in (jlOp (in units of the recoil velocity vq = Jcd/h) 
as a function of the phase shift in ([8]) for L = 16. The (to)-averaged velocity ([7]) obtained by the 
direct time propagation of the initial state up to 200T (dashed line) is compared to the asymptotic 
dependence given by the Floquet approach ^ (red solid line). Note the anti-symmetry behavior 
Vc{Q) = —Vc{@ + vr). The parameters are the same as in Fig. 2. 

direct time propagation. 

The resonance peaks can be associated with avoided crossings between two quasienergy 
levels (see Fig. These avoided crossings cause a strong velocity enhancement if one 

of the interacting, and transporting eigenstate overlaps significantly with an initial, non- 
transporting state of the motor. Note also that a very narrow avoided crossing requires a 
very large evolution time to become resolved, i.e., tobs ~ ^/|ea — ^fsl, see in Ref— . Clearly, 
our chosen evolution time of 200T is typically not large enough to clearly resolve the distinct 
resonances depicted in Fig. [21 

With Fig. m we show an example of the quasienergy spectrum, e„(0), and the carrier and 
starter velocities, Vc (solid line) and Vs (dashed line) respectively, for the coupling constant 
W/J = 0.01, and L = 4. Here we remark that a small coupling constant W yields the 
avoided crossings smaller while making the resonances sharper. This rather small system 
size, however, provides already a reasonable number of quasienergies for our elucidation. 
For every velocity resonance depicted in Fig. H] (bottom), one finds an avoided crossing in 
Fig. m (top). The fine avoided crossing structure induces an accompanying fine structure of 
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FIG. 4: Resonances of the motor velocity. The quasienergy spectrum (top) and asymptotic 
carrier /starter velocity (solid/dashed line) in units of vq = Jcd/fi as a function of 6 (bottom), 
obtained numerically for fiuj = 0.05 x J^,, Ai = 1, A2 = 0.5, W = 0.01 x J^, Jg = Jc, L = 4. 
Presence of the avoided crossings (emphasized by the circles and the rectangle frame) manifests in 
the resonant velocity peaks at the corresponding values of the phase 9. The insets (top) zoom into 
the region of one of the avoided crossing: the effect of a tiny anticrossing inside the larger one is 
clearly resolved in the sharp resonance on top of the wider one (bottom). 

corresponding velocity resonances, note the two insets in Fig. |H 

So far, we mainly focused on the motor velocity as given by the carrier subsystem velocity. 
Let us here briefly also comment on the starter dynamics, e.g., a possibly non-zero starter 
velocity. We found that the averaged starter velocity Vs indeed sensitively depends on the 
system parameters: It can either be very small compared to the carrier velocity (Fig. HI 
bottom, dashed line) or also larger than Vc- In short, the starter can move co-directionally 
or contra-directionally to the carrier motion. 

A robustness of Vcito) against a variation of to characterizes the quality of our quantum 
motor, as it has been mentioned in the abstract. Looking at the quasimomentum distribution 
(similar to that of on Fig. [2] (a,b)), it becomes evident that weekly depends on time. 
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FIG. 5: Dispersion of the motor velocity (jlip vs. the number of lattice sites L. The dependence 
of the asymptotic direct current on to vanishes with increase of the lattice size. The insets show 
the running average for the carrier velocity. For both insets two realizations of to, which give the 
maximal and minimal values of the asymptotic velocity are depicted (solid lines), with remaining 
realizations for different to varying in between these lines. The dependence averaged over 20 
realizations is depicted in both insets by a dashed line. Note that for L = 64 the range bounded 
by two curves practically shrinks to the single thick line. Here 6 = tt/2, and the other parameters 
are the same as in Fig. [2j 

therefore we could expect that the overlap with the initial state, c„(to), according to Q), and 
as a result i'c(^o) has a week dependence on to- To provide a more quantitative argument, 
we calculate the dispersion of the asymptotic motor velocity with respect to the switch on 
time to, 



Here we use the direct numerical time-evolution over the sufficiently long time, t = SOOT, 
to approach the asymptotic value ([7]), see insets in Fig. [5]). We found that the dispersion 
becomes increasingly negligible with increasing size L of the lattice. Starting out from 
L ^ 16, the motor gains practically the same asymptotic velocity independently on switch 
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FIG. 6: (color online) Average motor velocity vs. the load Asit) = UJBt. (a) The range of velocity 
values for the Floquet eigenbasis. (b) The average motor velocity for the initial condition with the 
localized carrier and the delocalized starter. We set ujb = uj-q/r, with r = 10, and vary the integer 
q. The parameters are W = 0.2Jc, Js = Jc = J, fiuj = O.lJc, Ai = 0.5, A2 = 0.25, = tt/2, and 
L = 4. 



on time to- 

This effect is due to the presence of the starter: The carrier velocity is obtained as the 
trace over the part of the total system Hilbert space, \lc)®\h), associated with starter degrees 
of freedom, \ls). The starter dynamics mimics a finite "heat bath" for the carrier dynamics 
whose effectiveness increases with both, the (i) the dimension of the starter subspace, i.e. 
the size L, and (ii) the strength of the interaction W. 



VI. LOAD CHARACTERISTICS 



Thus far, the analysis of the system in Eqs.([T]-[31 E]) has been performed in a idle-running 
mode with no load applied. In order to qualify for a genuine motor device, the engine 
must be able to operate under an applied load. The load is introduced in the form of an 
additional bias Asit) = ujst, being added to the vector potential A{t). All the information 
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about transport properties can be extracted by using again the Floquet formalism, provided 
that the ac-driving and the Bloch frequencies are mutually in resonance^^; i.e. we have the 
condition q ■ uo = r ■ ujb obeyed, where r and q are co-prime integers. Figure depicts the 
range of velocity values for Floquet eigenbasis, U„, see part (a), and the dependence of the 
resulting asymptotic motor speed, Vc-, for different bias values, part (b). There occur two 
remarkable features. First, the spectrum of velocities is symmetric around uob = 0. This 
follows because of the specific choice of the phase shift at O = 7r/2. Second, while some 
regimes provide a transport velocity along the bias, others correspond to the up-hill motion, 
against the bias. Therefore, a stationary transport in either direction is feasible. The load 
characteristics exhibits a discontinuous, fractal structure. In distinct contrast to the classical 
case^"^, it cannot be approximated by a smooth curve. This is a direct consequence of the 
above mentioned resonance condition. 



VII. EXPERIMENTAL REALIZATIONS 



We next discuss the parameter range where the carrier atom generates a tangible dc 
current. In the fast-driving limit, uj ^ iy/(/;.(|74i| + |v42|)), the Floquet states adiabatically 
follow the instantaneous eigenstates of the total Hamiltonian in presence of ([T]) zero particle- 
interaction i.e., W = Thus, the resulting, fast driven dc current approaches zero. In the 
slow-driving limit, u <^ iy/(^(|y4i| + |v42|)), the Floquet states adiabatically follow the 
associated instantaneous levels of the static Hamiltonian ([1]) with no dynamical symmetry- 
breaking field acting; thus the slow-driven dc current vanishes as well. Thus the maximum 
dc current region is located in the intermediate region, oj ~ + 1^21))- 

For an experimental realization of this quantum atom motor the following feature should 
be respected: Because in the tight-binding approximation the maximal amplitude of the 
tunneling is limited from above, Jc ^ Jmax = CIS-Eq, for, e.g. a ^Li "carrier", a lattice 
constant d ~ 10 fim, and with hu = O.lJc as used in the present calculations, the driving 
frequency u should be less than 2Hz. Then, the time required to launch the motor; i.e., for 
it to approach the asymptotic velocity, which is ~ 0.05^0 ~ 30fim/s ~ 3 sites/second, is 
around a minute. Further focusing of the laser beam can decrease the lattice constant d, 
thereby decreasing the launch time to experimentally accessible coherence times around 10 
seconds^ and increasing the asymptotic velocity. 
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VIII. CONCLUSIONS 



We elaborated in greater detail on the working principles of the quantum electric ac-motor 
that is made of two ultracold interacting atoms only^^: a carrier and a starter, moving in a 
optical bracelet potential. The central result of this study is an evident directed coherent 
carrier motion that is induced by the starting mechanism. The emerging motor velocity can 
suitably be controlled by means of the symmetry breaking, time-dependent, bichromatic, 
external magnetic flux. Importantly, for zero-momentum initial conditions the asymptotic 
carrier velocity loses its dependence on the switch-on time to of the ac-drive upon increasing 
the bracelet size L. 

An extension of our motor setup to several interacting carries or starters (i.e., multiple 
rotor motors or finite bosonic "heat baths") presents an intriguing challenge. A particular 
interesting objective to pursue is the problem of whether the motor velocity can be optimally 
tuned with the number of participating atoms? 

Finally, an exciting perspective is to physically couple our quantum motor to a nano- 
mechanical resonator^'^. Such a hybrid system can be used to power quantum mechanically 
such a classical object. 

This work was supported by the DFG through grant HA1517/31-1 and by the German 
Excellence Initiative "Nanosystems Initiative Munich (NIM)". 
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